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Abstract 



The yrast spectra (i.e. the lowest states for a given total angular momen- 
tum) of quantum dots in strong magnetic fields, are studied in terms of exact 
numerical diagonalization and analytic trial wave functions. We argue that 
certain features (cusps) in the many-body spectrum can be understood in 
terms of particle localization due to the strong field. A new class of trial 
wavefunctions supports the picture of the electrons being localized in Wigner 
molecule-like states consisting of consecutive rings of electrons, with low-lying 
excitations corresponding to rigid rotation of the outer ring of electrons. The 
geometry of the Wigner molecule is independent of interparticle interactions 
and the statistics of the particles. 

PACS: 67.40.Db, 73.21.La 
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I. INTRODUCTION 



In a sufficiently strong magnetic field, the spin degree of freedom of electrons in quantum 
dots is frozen out, and they may be treated as spinless. Furthermore, if the electron density is 
low, the strong field will cause localization of the electrons in a Wigner molecule (a finite size 
realization of a Wigner crystal). For a circular dot, this localization happens gradually as the 
magnetic field increases and is ffist seen as a separation of a ring of electrons from the 
so-called maximum density droplet (MDD) , the finite-size analog of the u = 1 quantum 
Hall 1^ state with angular momentum N{N—l)/2. By increasing the field further, eventually 
all the electrons localize as in a classical one-component plasma. This localization takes 
place in the internal coordinate system of the dot and may be revealed by using a rotating 
frame or studying pair correlation functions |^. 

In this paper we study the interplay between particle localization in a finite 2D electron 
system in a harmonic confinement and the structure of its yrast spectrum, i.e. its lowest- 
lying many-body states as function of the total angular momentum L. Several classes of 
trial many-body wavefunctions are available to describe the yrast states of such a system. In 
particular, electron states with angular momenta equal to odd multiples of N{N — l)/2 can 
be approximately described by the Laughlin wavefunction for the fractional Hall effect [§. 
Another class of wavefunctions that has proven successful in describing yrast states away 
from these special L- values P,|10| is due to Jain [pJ]] . 



Here, we introduce an additional type of trial wavefunctions which, although perhaps 
not as accurate as the Jain construction, in our opinion give a more intuitive picture of the 
structure of the ground state of the dot as one goes from a Laughlin state, L = qN{N — l)/2 
to the quasihole- or single vortex state L = qN{N — l)/2 + N. These wavefunctions are 
similar to those recently introduced by Bertsch and Papenbrock (BP) [|12| to describe the 
ground state of a weakly repulsive, harmonically trapped, rotating 3D Bose condensate for 
2 < L < N. The close analogy between this system and a fermionic quantum dot in the 
quantum Hall regime has been discussed in a number of recent papers [|T3|-0; the essential 
observation is that the weak interaction limit allows for a two-dimensional description of 
the boson system in terms of lowest Landau level (LLL) wavefunctions which, just as for 
the quantum Hall effect, may contain powers of the Laughlin- Jastrow factor Yl^^j{zi — zj), 
with Zi the complex coordinate of the ith particle. The only difference is the symmetry of 
the many-body wavefunction. This means that it is possible to 'borrow' techniques familiar 
from quantum Hall physics to describe the rotating Bose condensates and vice versa. 

The outline of the paper is as follows. In section 2 we present exact numerical calculations 
for our test case, the six-electron dot. Trial wavefunctions are introduced and compared to 
the exact results in section 3. Finally, we discuss the many-body spectra in terms of electron 
localization, using both our numerical and analytical results. 



II. EXACT CI CALCULATIONS AND THE YRAST SPECTRA 

Independent electrons in a strong magnetic field occupy Fock-Darwin states |lT7| , p!8| with 
the effective confinement potential determined by the frequency = \JojI + oj^/A, where 
ujq is the external confinement and uj^ the cyclotron frequency eB/m. In the absence of 
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interactions, the ground- and low-lying excited states (for high angular momentum values) 
are entirely in the lowest Landau level (LLL). 

In order to solve the many-body problem in the presence of Coulomb interactions, we 
use the configuration interaction technique (CI); it has the advantage that the many-body 
energy spectrum, not only the ground state, can be computed and that in the case of small 
particle numbers the solution is essentially exact. We choose a six electron quantum dot as 
our test case. This is still so small that in the presence of a magnetic field (the electron 
system being polarized) an essentially exact result can be obtained for several low-energy 
many-body states. The six electron system is the smallest size where the classical Wigner 



molecule has two nearly degenerate isomers [^], and as we shall see, the interplay of these 
isomers leads to interesting features in the many-body spectrum, as first pointed out by 
Maksym 0|. 

In the CI computations we use the Darwin-Fock states (with frequency Ucs) as our 
single particle basis. For the smallest fields, where the ground state is the MDD, it was 
essential to include single particle states from several Landau levels. The Fock states included 
were chosen with the help of a cut-off procedure: The matrix dimension was restricted to 
about 30000, which amounts to including about 50 Fock states. In the high-field case, 
corresponding to filling factor 1/3, the single particle basis was restricted to the LLL, and 
a complete CI was used with single particle states up to £ = 20. (Note that the g = 3 
Laughlin state, to be discussed below, only includes single particle states up to ^ = 15 for 
six electrons). 

We calculated the many-body spectra for three different values of the magnetic field 
{b = 0.9, 1.2 and 2.5 in atomic units, corresponding approximately to filling factors 1, 2/3 
and 1/3, respectively). The strength of the external confinement was ujq = 0.34435 a.u. 
(corresponding to a 2D electron density at the equilibrium density, = 1.51 a.u.). The 
many-body spectra are shown in Figs. |I] and ^ For each field strength we show the original 
spectrum (Fig. 0) and a spectrum displaying the interaction energy, where for each angular 
momentum value the center of mass contribution Lhues has been substracted (Fig. 

The results show that the overall structure of the spectrum is independent of the magnetic 
field, the role of which is mainly to tilt the spectrum so that the minimum energy will be at 
a finite angular momentum pOf. The spectrum starting from Li = 15 is nearly identical for 
the field strengths 0.9 and 1.2 and very similar to the spectrum corresponding to the strong 
field 6 = 2.5 after the angular momentum L3 = 45. 



III. TRIAL WAVEFUNCTIONS 

In order to get some understanding of the microscopic structure of the yrast states, it is 
useful to study them in terms of trial wavefunctions. 

At angular momenta Lg = qN{N — l)/2 (with q an odd integer), corresponding to filling 
fraction u = 1/q, an obvious guess for the ground state is the Laughlin wavefunction |^, 

N 

V'L(g) = n(^^-^.)'e-^^^'''^''^ (1) 

i<j 

where Zk = JmiUcs/q{xk + iyk) are the complex particle coordinates. 
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For q = 1 this wavefunction is a Slater determinant of LLL single particle states, rjg oc 
exp{—qzl/4), corresponding to angular momenta from i = to i = N — 1. It has been 
shown that the wavefunction (|I]) describes very well the exact many-body state when the 
angular momentum corresponds to g = 3 or larger g-values, both in a spherical geometry 
(which mimicks a homogeneous state) El]] and for a disk geometry with "open" boundary 



conditions ||8|,|22[|. In the case of g = 1 the Laughlin wavefunction reduces to the (restricted) 
Hartree-Fock approximation which, as we shall see, is not a quantitatively exact description 
to the state for = 6, as the latter contains significant contributions from higher Landau 
levels. 

For angular momenta above those described by the Laughlin wavefunction, L = Lq + M, 
ansatz wavefunctions can be constructed by multiplying the Laughlin- Jastrow factor by a 
fully symmetric polynomial in the Zk of order M. In this way, the good correlations of the 
Jastrow factor are still built in. For M = 1, the only way of doing this is multiplying hjj^^k- 
This is just a center-of-mass excitation of the Laughlin state. (In a harmonic confinement 
the center of mass motion separates from the internal motion.) For 2 < M < N we propose 

M 

i^LBp{Lg + M)=S l[{z, - zo) i^Liq) (2) 

i 

where the operator S symmetrizes the product, and zq = J^^k/N denotes the center of mass 
coordinate. 

The first, symmetric part of this wavefunction was introduced by Bertsch and Papen- 



brock ||12[ to describe the ground state of a weakly interacting, rotating Bose condensate for 
angular momentum up to the "unit vortex" state, i.e. L < N. In their case of a delta func- 
tion interaction, the Bertsch-Papenbrock wavefunction has been proven to be exact p3|-|25 



Here, we combine it with the Laughlin- Jastrow factor to construct a (fully antisymmetric) 
wavefunction that interpolates between the Laughlin state L = Lg and the "quasihole" state 
L = Lq + N . Note that, indeed, for M = Eq-O) just takes the form of the wavefunction 
proposed by Laughlin 0, for a quasihole centered at zq. 

We have tested the applicability of the Laughlin-Bertsch-Papenbrock (LBP) wavefunc- 
tion (^ in the present case of six Coulomb-interacting electrons in a harmonic potential. 

Table 1 gives the calculated overlaps of this wavefunction with the exact wavefunction 
starting from Li = 15 and from L3 = 45. In the first case we determined only the values 
for Li -|- 5 and Li + 6 which are related to the localization of electrons (see next section). 
The LBP wavefunction seems to be a fairly good approximation also for the fermion case, 
especially for g > 3; the reason that the overlaps for L = 20, 21 are relatively poor, is 
that the Jastrow factor itself is not a very good approximation to the exact solution at 
L = 15. As discussed previously, this is due to contributions from higher Landau levels in 
the exact wavefunction. If we perform an exact diagonalization for q = 1 but restricted to 
the LLL, these latter overlaps will be much better (1.00, 0.96 and 0.99 for L = 15, 20 and 
21, respectively). 

In Table 1 we also give the most important single particle configuration of the CI calcu- 
lation and its amplitude. The MDD is strongly dominated by only one configuration while 
the L3 case is much more dispersed. We have argued that the Laughlin-Bertsch-Papenbrock 
ansatz cannot be accurate after the MDD; nevertheless, the most important configuration 
is consistent with the qualitative picture given by the LBP ansatz: For each additional 
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quantum of angular momentum, one particle is pushed outwards from the center of the dot 



There is another class of trial wavefunctions |Tl|, based on the idea of mapping strongly 
interacting electrons (or bosons) onto weakly interacting composite fermions by attaching 
an even (odd) number of flux quanta to each particle. These so-called Jain states have 
proven very successful in describing yrast states in quantum dots |jl0|,^,|23 as well as weakly 
interacting rotating Bose condensates [|T3|,|TB|. Especially, the method explains the main 
features of the yrast spectra already in the noninteracting illustrated in Fig. H, and 

(in the interacting case) gives excellent overlaps with exact wavefunctions. However, the 
Jain construction involves technical steps (formally constructing states in higher Landau 
levels and subsequently projecting down to the LLL in a not uniquely defined way), which 
generally make it hard to visualize the qualitative features of the state. For this reason we 
choose to study the LBP ansatz which might not be as accurate as the Jain construction, 
but which in a very intuitive way "interpolates" between the Laughlin state Lq and the 
'quasihole' state Lq + N and in a transparent way links the many-body spectrum to Wigner 
crystallization. 



IV. LOCALIZATION OF ELECTRONS AND EDGE RECONSTRUCTION 

The original spectra of Fig. |l|show energy minima (cusps) at angular momenta 15, 20, 25, 
30 (6 = 0.9 and 1.2) and at 45, 50, 55, 60 (6 = 2.5). This can be understood by localization 
of electrons in a five- fold ring (with one atom at the center). A rigid rotation of a ring of N 
electrons can only have angular momenta L = nN if is odd and L = [n + 1/2)N if is 
even (n is an integer) [^. Other angular momentum values can only be a result of center 
of mass excitations or as a combination of the rigid rotation with vibrational modes. The 
analysis of rotation- vibration spectra was found to be successful also in studying localization 



of non-polarized electrons in quantum rings |29 



The spectra also show markedly low energies when L is 21, 27, 51, 57, corresponding to a 
six-fold symmetry. Classically the six-electron system has also a stable geometry of a six-fold 
ring |]19[, which is only slightly higher in energy than the five-fold ring. In the special points 



L = 15 and L = 45 the quantum state must be a combination of these two structures. If 
the angular momentum is increased by 5, the solution collapses into the five-fold symmetry, 
and if it is increased by 6 the solution collapses into the six-fold ring. 

Cusps of the type seen here, have previously been noticed for a few-electron system and 
discussed in terms of particle localization |20| and symmetry |3yj31| as well as composite 



fermions |]rU|. Here we will first use the exact solutions for the six electron dot to test our 



qualitatively simple picture of localized electrons. 

Figures |^ and ^ show the radial densities of the CI computations for different values of 
the angular momentum at the yrast line. Clearly the densities for L = 20, 25, 45 and 50 
are consistent with the five-fold ring, having one particle in the center, while the results for 
L = 21, 27, 51 have only one ring of electrons. 

It is illustrative to study the electron localization also in terms of analytic trial wave- 
functions. The W-particle correlation function' is simply defined as the square of the wave- 
function 
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p{zi, Z2, ■ ■ ■ , zn) = \i^{zi, Z2,--- , Zn)\'^- (3) 

Laughlin [§] showed that in the case of Eq. finding the maximum of the correlation func- 
tion is equivalent to finding the equilibrium electron positions in a classical one-dimensional 
plasma, where the electrons interact via a repulsive logarithmic potential and are confined by 
a harmonic trap. In accordance with the above discussion, for six particles and L = 15 or 45, 
there are two local minima corresponding to the five- and six-fold geometries, respectively 
(the former being the lowest). 

Moreover, Laughlin showed that the optimal geometry is independent of the exponent 
q of the wavefunction. Scaling the complex coordinates with q as we did in Eq. (|T|), also 
the distances between the localized electrons are independent of q. The Laughlin wave 
function for even g- values describes well rotating Bose condensates [ll5| , p!6| . Consequently, 



at large angular momenta the statistics of the particles do not play any role in the internal 
correlation of the particles. In fact, this picture is also fairly independent of the actual 
interaction potential, since the Laughlin state is an excellent approximation for long-range 
Coulomb interactions as well as zero-range interactions 0. 

The localization becomes stronger when q increases, making the Laughlin state gradually 
more similar to a crystal 0. This can be measured by the sharpness of the maximum of p. 
It is easy to show that around a maximum correlation 



\dxl ' / M^i' ^2,- ■■,^iv) 1,^^,0 = -2gp(2i, 2:2, ■■■,2;v), (4) 

where are the coordinates of the maximum. This result explains why the density of the 
MDD is nearly constant, while the density corresponding to g = 3 has sharper maxima 
consistent with a configuration of electrons localized in a five-fold geometry. 

Also, the exact many-body spectra of Fig. |l| show clearly that the five-fold ring becomes 
more dominating after L3 than after Li, as it drops below the six-fold one for all L > L3. 

Our analytical ansatz of Eq. (^) allows for the study of localization of electrons at any 
angular momentum after any Lg. Figure ^ shows the positions of electrons at the correlation 
maxima of Eq. (|]) for L = 45, 47, 48, 49, 50, and 51 {L = 46 is an uninteresting center 
of mass excitation). The five- fold ring at 50 and the six- fold ring at 51 are reproduced as 
anticipated. In general, the outer ring consists of as many electrons as is the order of the 
polynomial multiplying the Laughlin wavefunction. 

Note that the geometries are independent of g; for example they would be the same for 
bosons at L values ranging from 30 to 36. (Or, in fact, for non-integer q, which corresponds 
to anyons, i.e. particles obeying fractional statistics, in the LLL p2| ) 

For systems containing a large number of electrons, it is not possible to perform exact 
calculations of the yrast spectra. However, in the region where the particles localize, simple 
arguments can be used to predict properties of the yrast spectra. The localization in the 
Laughlin state corresponds to the classical plasma with logarithmic interactions. The ground 
state positions of the electrons are very similar to those obtained with 1/r interactions, 
which have been analyzed by Schweigert and Peeters [^. In finite systems the resulting 
Wigner molecules are characterized by rings of electrons (as seen already in the case of six 
electrons). Figure |^ shows as an example the 18 electron geometry. If the Laughlin state 
is again multiplied by a polynomial as in Eq. (0), it follows that especially favourable are 
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cases where the order of the polynomial equals the number of electrons in the outermost 
ring. This case is demonstrated in Fig. ^ which shows that when the angular momentum 
is increased by 11, i.e. the number of electrons in the outer ring, the only effect is the 
expansion of this outer ring. 

The above analysis suggests that starting from the Laughlin state, the next low-energy 
state (cusp) is obtained at angular momentum Lq+n where n is the number of electrons in the 
outermost ring of the corresponding classical electron system. In the case of the maximum 
density droplet this corresponds to the expansion of the so-called Chamon-Wen Q edge 
from the MDD. Indeed the number of electrons in the Chamon-Wen edge as calculated by 
Reimann et al. for = 20 is in agreement with the number of classical electrons in the 
outermost ring. Increasing the angular momentum further by the number of electrons in 
the next ring causes also that ring to expand, and again an energy minimum is expected. 

The present scheme thus provides microscopic trial wavefunctions describing the edge 
reconstruction of small quantum Hall droplets. For the special case of the MDD (g = 1) a 
similar trial wavefunction was recently proposed by Goldmann and Renn 0]. 

The expansion of the Chamon-Wen ring also increases the electron localization. This can 
be demonstrated by studying the decrease of the A^-particle correlation when an electron 
is moved out from the position of the maximum correlation (keeping all other electrons 
in place). Figure shows the decrease of the correlation as a function of the electron 
displacement (angular average) for the six electron dot for the MDD (L = 15), for the 
Chamon-Wen state {L = 20) and for the g = 3 Laughlin state {L = 45) calculated using the 
LBP ansatz Eq.(0). It is clear that in the Chamon-Wen edge state the electron localization 
is increased as compared to the MDD. The difference between the L = 45 and L = 15 states 
in Fig. ^is in agreement with Eq. (^. 

Let us finally comment on the case L = Lg + N. Our ansatz (|^) then has the same form 
as that proposed by Laughlin to describe the fractionally charged quasihole excitations 
of large quantum Hall systems. In these systems, the quasihole is the lowest-lying bulk 
excitation. The finite size of the quantum dot, with no clear distinction between bulk and 
edge, allows also for low energy excitations where one surface layer (of localized electrons) at 
a time is excited. As we have argued, the classical interpretation of these excitatitons is that 
the outer ring is rigidly rotating. The lowest energy excitation then corresponds to angular 
momentum increase which equals the number of electrons in the outermost ring. This is 
already seen in the six electron case where the low energy excitations after the L3 = 45 state 
have L = 50, 55 etc in accordance with the 5 electrons in the ring around the center atom. 



V. CONCLUSIONS 

We have demonstrated that the localization of electrons in a strong magnetic field can 
be seen in the structure of the many-body spectrum. The many-body yrast states following 
each Laughlin state can be approximated with simple trial wavefunctions similar to those 
originally suggested by Bertsch and Papenbrock |T2| for rotating Bose condensates. 



In a quantum dot the internal structure of the Laughlin state corresponds to localization 
of electrons in consecutive rings. The Laughlin-Bertsch-Papenbrock wavefunction illustrates 
that the Chamon-Wen edge reconstruction can be seen as an expansion of the outermost 
ring. Since the Laughlin and LBP wavefunctions are good approximations for fermions as 
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well as bosons and both for long- and short range interactions, it seems that the internal 
geometry of the localized particles in a high rotational state depends only on the number of 
particles in the harmonic confinement. 
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FIGURES 



FIG. 1. Calculated many-body states as a function of angular momentum for three different 
values of the external magnetic field. The lowest energy states corresponding to a five-fold (six-fold) 
symmetry are connected with a dashed (dotted) line. 

FIG. 2. Many-body states of Fig. |^ but with the center of mass excitation energy LhtOeS 
substracted. The structure of the yrast line predicted by the Jain construction is shown as dashed 
lines. 

FIG. 3. Radial density distributions of the exact many-body states in a six electron quantum 
dot at a magnetic field of 6 = 2.5 a.u. L is the total angular momentum. 

FIG. 4. Radial density distributions of the exact many-body states in a six electron quantum 
dot at a magnetic field of 6 = 1.2 a.u. L is the total angular momentum. 

FIG. 5. Geometry of the 'Wigner molecule' of six electrons for different values of the angular 
momentum L, calculated from the trial wavefunction Eq.(|2|). 

FIG. 6. Positions of electrons in an 18-electron dot at the maximum of the square of the 
wavefunction (^) for L = qN{N — l)/2 = 153g (filled circles) and L = 153q + 11, 11 being the 
number of electrons in the outermost ring. 

FIG. 7. The decrease of A^-particle correlation when one electron is moved out from the 
maximum position. The distance is given as percentage of the nearest neighbour distance. Solid 
line L = 15, dashed line L = 20, dotted line L = 45. The dot has six electrons. 
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TABLES 



TABLE L Absolute values of the overlaps between the exact wavefunction and the approximate 
wavefunction, Eq. (|2[), for different values of the angular momentum L. MIC is the most important 
single particle configuration as occupation numbers of the LLL states starting from i = 0. The 
last column gives the amplitude of the MIC. 
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